Theory of the leak-rate of seals 
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Seals are extremely useful devices to prevent fluid leakage. However, the exact mechanism of 
roughness induced leakage is not well understood. We present a theory of the leak-rate of seals, which 
is based on percolation theory and a recently developed contact mechanics theory. We study both 
static and dynamics seals. We present molecular dynamics results which show that when two elastic 
solids with randomly rough surfaces are squeezed together, as a function of increasing magnification 
or decreasing squeezing pressure, a non-contact channel will percolate when the (relative) projected 
contact area, A/Aq, is of order 0.4, in accordance with percolation theory. We suggest a simple 
experiment which can be used to test the theory. 



1. Introduction 

A seal is a device for closing a gap or making a joint 
fluid tight [ij. Seals play a crucial role in many modern 
engineering devices, and the failure of seals may result in 
catastrophic events, such as the Challenger disaster. In 
spite of its apparent simplicity, it is still not possible to 
predict theoretically the leak-rate and (for dynamic seals) 
the friction forces[2j for seals. The main problem is the 
influence of surface roughness on the contact mechanics 
at the seal-substrate interface. Most surfaces of engi- 
neering interest have surface roughness on a wide range 
of length scales [1], e.g, from cm to nm, which will influ- 
ence the leak rate and friction of seals, and accounting for 
the whole range of surface roughness is impossible using 
standard numerical methods, such as the Finite Element 
Method. 

In this paper we will analyze the role of surface rough- 
ness on seals. We will use a recently developed contact 
mechanics theoryjl, [1,0,0, [1,0] to calculate the leak-rate 
of static seals. We assume that purely elastic deforma- 
tion occurs in the solids, which is the case for rubber 
seals. For metal seals, strong plastic deformation often 
occurs in the contact region. 




FIG. 1: Rubber seal (schematic). The liquid on the left-hand- 
side is under the hydrostatic pressure Pa and the liquid to the 
right under the pressure Pb (usually, Pb is the atmospheric 
pressure). The pressure difference AP = Pa — Pb results in 
liquid flow at the interface between the rubber seal and the 
rough substrate surface. The volume of liquid flow per unit 
time is denoted by Q, and depends on the squeezing pressure 
Po acting on the rubber seal. 



The theory developed below is based on studying the 
interface between the rubber and the hard countersur- 
face (usually a metal) at different magnifications C- At 
low magnification the surfaces appears flat and the con- 
tact between them appears to be complete (i.e., no leak 
channels can be observed). However, when we increase 
the magnification we observe surface roughness at the 
interface, and, in general, non-contact regions. As the 
magnification increases, we will observe more and more 
(short-wavelength) roughness, and the (apparent) con- 
tact area A{C,) between the solids will decrease. At high 
enough magnification, for C = a non-contact (percola- 
tion) channel will appear, through which fluid will flow, 
from the high pressure side (pressure Pa) to the low pres- 
sure side (pressure Pb), see Fig. [T] We denote the most 
narrow passage between the two surfaces along the perco- 
lation path as the critical constriction. When the magni- 
fication increases further, more percolation channels will 
be observed, but these channels will have more narrow 
constrictions than those for the first channel which ap- 
pears at the percolation threshold (C = Cc)- 

The picture described above for the leakage of seals has 
already been presented by one of the present authors [1, 
\ld^ . However, recent development in contact mechanics 
now allows us to present a more accurate analysis of the 
leaka ge p rocess. In this paper we extend the theory of 
Ref. [lO|, and present numerical results for the size of 
the critical constriction and for the leak-rate. 

In Sec. 2 we describe the basic picture used to calculate 
the leak-rate of static seals. In Sec. 3 we present numer- 
ical results for the size of the critical constriction and for 
the leak-rate. In Sec. 4 we present Molecular Dynam- 
ics results which illustrate how the contact between the 
two solids changes as the magnification C increases. We 
find that the percolation channel is formed when C = Cc, 
where A(Cr)Mn ~ 0.4, in accordance with percolation 
theory In Sec. 5 we improve the theoretical picture 
of how to understand static seals. In Sec. 6 we compare 
the theory with experimental data. In Sec. 7 we present 
some comments related to the non-uniform seal pressure 
distribution, the role of adhesion and rubber viscoelas- 
ticity. In Sec. 8 we study dynamical (linear reciprocal 



FIG. 2: The rubber-countersurface apparent contact area is 
rectangular Ly. We "divide" it into A'^ = Ly/L^ square 

areas with side L = and area Ao = . 




FIG. 3: An rubber block (dotted area) in adhesive contact 
with a hard rough substrate (dashed area). The substrate 
has roughness on many different length scales and the rubber 
makes partial contact with the substrate on all length scales. 
When a contact area is studied at low magnification it appears 
as if complete contact occurs, but when the magnification is 
increased it is observed that in reality only partial contact 
occurs. 

motion) seals at low sliding velocities. In Sec. 9 we sug- 
gest a simple experiment to test the theory. Sec. 10 
contains the summary and the conclusion. 

2. Theory 

We first briefly review the basic picture on which our 
calculations of the leak-rate is based Assume that 
the nominal contact region between the rubber and the 
hard countersurface is rectangular with area Lx ^ sqb 
Fig. O We assume that the high pressure fluid region 
is for a; < and the low pressure region for x > L^- 
We now divide the contact region into squares with the 
side Lx — L and the area = (this assumes that 
N = Ly/Lx is an integer, but this restriction does not 
affect the final result). Now, let us study the contact 
between the two solids within one of the squares as we 
change the magnification We define Q ^ L/X, where 
A is the resolution. We study how the apparent contact 
area (projected on the xy-plane), ^(C), between the two 
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FIG. 4: The contact region at different magnifications 
(schematic). Note that at the point where the non-contact 
area (white area) percolate A{C,c) ~ 0.4Ao, while there ap- 
pear to be complete contact between the surfaces at the lowest 
magnification C = 1- ^(1) = ^o. 

solids depends on the magnification C,. At the lowest mag- 
nification we cannot observe any surface roughness, and 
the contact between the solids appears to be complete 
I.e., A{\) = Aq. As we increase the magnification we will 
observe some interfacial roughness, and the (apparent) 
contact area will decrease, see Figs. |3] and ID At high 
enough magnification, say C = Cc, a percolating path of 
non-contact area will be observed for the first time, see 
Fig. m The most narrow constriction along the perco- 
lation path will have the lateral size Ac = i/Cc aud the 
surface separation at this point is denoted by Uc = wiCCc) 
and is given by a recently developed contact mechanics 
theory (see below). As we continue to increase the mag- 
nification we will find more percolating channels between 
the surfaces, but these will have more narrow constric- 
tions than the first channel which appears &i C, = and 
for the moment we will neglect the contribution to the 
leak-rate from these channels (see also Sec. 5). Thus, in 
this section we will assume that the leak-rate is deter- 
mined by the critical constriction. 

A first rough estimate of the leak-rate is obtained by 
assuming that all the leakage occurs through the critical 
percolation channel, and that the whole pressure drop 
AP = Pa — Ph (where Pa and Pb is the pressure to the left 
and right of the seal) occurs over the critical constriction 
[of width and length Ac ~ i/Cc and height Uc — ui(Cc)]- 
Thus for an incompressible Newtonian fiuid, the volume- 
flow per unit time through the critical constriction will 
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magnification C, 



Q = MAP, 



where 
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where 77 is the fluid viscosity. In deriving (1) we have 
assumed laminar flow and that Uc << Ac, which is al- 
ways satisfied in practice. Here we have introduced a 
factor a which depends on the exact shape of the critical 
constriction, but which is expected to be of order unity. 
Since there are N = Ly/L^ square areas in the rubber- 
countersurface (apparent) contact area, we get the total 
leak-rate 
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FIG. 5: An asperity contact region observed at the magnifica- 
tion (,. It appears that complete contact occur in the asperity 
contact region, but upon increasing the magnification it is ob- 
served that the solids are separated by the average distance 
u(C). 



To complete the theory we must calculate the separa- 
tion Uc = Mi(Cc) of the surfaces at the critical constric- 
tion. We first determine the critical magnification by 
assuming that the apparent relative contact area at this 
point is given by site percolation theory. Thus, the rel- 
ative contact area A{()/A() « 1 — Pc, where Pc is the so 
called site percolation thresholdfll]. For an infinite-sized 
systems Pc « 0.696 for a hexagonal lattice and 0.593 for a 
square lattice [ll|. For finite sized systems the percolation 
will, on the average, occur for (slightly) smaller values of 
p, and fluctuations in the percolation threshold will occur 
between different realization of the same physical system. 
We will address this problem again later (see Sec. 4) but 
for now we take pc w 0.6 so that ^(Cc)/^o ~ 0.4 will 
determine the critical magnification — (c- 

The (apparent) relative contact area A{C,)/Aq at the 
magnification C can be obtained using the contact me- 
chanics formalism developed elsewhere [3, S 0, B 
where the system is studied at different magnifications 
C, see Fig. [3l We havcjj^ ^] 
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where the surface roughness power spectrum 
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where (...) stands for ensemble average. Here E and v 
are the Young's elastic modulus and the Poisson ratio of 
the rubber. The height profile h{x) of the rough surface 



can be measured routinely today on all relevant length 
scales using optical and stylus experiments. 

We define ui(C) to be the (average) height separating 
the surfaces which appear to come into contact when the 
magnification decreases from C, to ( — A^, where AC is a 
small (infinitesimal) change in the magnification. ui(C) 
is a monotonically decreasing function of C, and can be 
calculated from the average interfacial separation u{C) 
and A{() using (see Ref. [9j) 

m{0=u{0+u'{OA{0/A'{C)- 

The quantity u{C) is the average separation between the 
surfaces in the apparent contact regions observed at the 
magnification C, see Fig. [5l It can be calculated fromQ 

S(C) - A r dq q'C{q)w{q) 



r°° 1 

X / dp' - [7 + 3(l-7)P'(g,y,C)] 
JpiC) P' 

where 7 « 0.4 and where 

p(C) = PoAo/A{0 

and 



q X -1/2 



wiq,0={7T dq' q'''C{q') 
The function P{q,p, () is given by 



2 /•s(9,C)p 

P{q,P,0 = —r dx e- 

Vi" Jo 

where s{q, () = w{q, Q/E* . 



We study the contact between the sohds at increas- 
ing magnification. In an apparent contact area observed 
at the magnification C, the substrate has the root mean 
square roughness ampUtude[3. \ld^ 

h'^UO ^ 27r r dqqCiq). (4) 
-'Cqo 

When we study the apparent contact area at increasing 
magnification, the contact pressure will increase and 
the surface roughness amplitude /irms(C) '^iH decrease. 
Thus, the average separation u(C), between the surfaces 
in the (apparent) contact regions observed at the magni- 
fication will decrease with increasing magnification. 

3. Numerical results 

We now present numerical results to illustrate the the- 
ory developed above. We assume a rubber block with a 
flat surface, squeezed by the nominal pressure Pq against 
a hard solid with a randomly rough surface which we as- 
sume to be self affine fractal. Thus the surface roughness 
power spectrum for qq < q < qi: 

C{q)^Coq-'^'+"^ 

where 

where go and qi are the long-distance and short-distance 
cut-off wavevectors, respectively. 

The rubber has the Young's modulus E ~ 10 MPa 
(as is typical for the low-frequency modulus of rubber 
used for seals) and Poisson ratio v = 0.5. The pressure 
difference in the fluid between the two sides of the seals is 
assumed to be AP = 0.01 MPa. The fluid is assumed to 
be an incompressible Newtonian fluid with the viscosity 
7] = 0.001 Ns/m^. We will study how the lateral size Ac 
and the height Uc of the critical constriction depends on 
the fractal dimension Df = 3 — H and on the root-mean- 
square roughness amplitude h^ms of the rough surface. 
We also present results for how the volume flow of fluid 
through the seals depends on D{ and /irms- The randomly 
rough surfaces have the cut-off wavevectors qo — 1.0 x 
10^ m~^ and qi = 7.8 x 10^ m~^, and we vary the applied 
squeezing pressure Pq from 0.05 MPa to 1 MPa. 

Let us first vary the rms roughness amplitude. In Fig. 
[6] we show the lateral size Ac = A(Cc) and in Fig. [7] the 
height (interfacial separation) Uc of the critical constric- 
tion, as a function of the applied normal (or squeezing) 
pressure Pq ■ Results are shown for self affine fractal sur- 
faces with the Hurst exponent H = 0.8 (or fractal di- 
mension Df — 2.2), and for surfaces with the root-mean- 
square roughness (rms) 1,2,4 and 6 /xm. As expected, 
the size of the critical constriction increases when the 
roughness increases. In Fig. [8] we show the volume per 
unit time of fluid leaking through the seals as a func- 
tion of the applied normal (or squeezing) pressure Pq. 
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FIG. 6: The lateral size Ac ~ ^{(c) of the critical constriction 
of the percolation channel, as a function of the applied normal 
(or squeezing) pressure Pq. Results are shown for self affine 
fractal surfaces with the Hurst exponent H = 0.8 (or fractal 
dimension Dt — 2.2), and for surfaces with the root-mean- 
square roughness (rms) 1, 2, 4 and 6 ^m. 

Note the extremely strong decrease in Q with increasing 
squeezing pressure and also its strong dependence on the 
rms roughness amplitude. 

In Figs. [HI Uni and [11] we show the analogous re- 
sults when we vary the Hurst exponent H — 0.9, 0.8, 
0.7 and 0.6 for hrms = 2 /j,m. Note that when H de- 
creases for a fixed /irms, the short wavelength roughness 
increases while the long wavelength roughness is almost 
unchanged. 

In Fig. [12] we show the interfacial separation ui(Cc) 
and the rms-roughness ft.rms(Cc) in the critical constric- 
tion, as a function of the applied normal (or squeezing) 
pressure Pq. Results are shown for a self affine fractal sur- 
face with the Hurst exponent H = 0.8 (or fractal dimen- 
sion Df — 2.2), and with the root-mean-square roughness 
(rms) 6 ^m. Note that the difference between /irms(Cc) 
and Mi(Cc) is relatively small. We have found that this is 
the case also for the other parameters used in the study 
above. 

4. Molecular Dynamics results 

The multiscale molecular dynamics model has been 
described in Ref . [l3| . but we review it briefly here. In 
what follows we denote the lower solid as substrate, the 
upper solid as block. We are concerned with the con- 
tact between a randomly rough and rigid substrate, and 
an elastic block, without adhesion. We are interested in 
surfaces with random roughness with wave-vector com- 
ponents in the flnite range qi > q > qL (see Fig. rT3| . 
where qL = 27r/L, L is the lateral size of the system. 
In order to accurately study contact mechanics between 
elastic solids, it is necessary to consider a solid block 
which extends a distance '--^ L in the direction normal 
to the nominal contact area. This requires huge number 
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FIG. 7: The interfacial separation Uc = mi(Cc) at the critical 
constriction of the percolation channel, as a function of the 
applied normal (or squeezing) pressure Pq. Results are shown 
for self afiine fractal surfaces with the Hurst exponent H = 0.8 
(or fractal dimension Df — 2.2), and for surfaces with the 
root-mean-square roughness (rms) 1,2,4 and 6 fim. 
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FIG. 8: The volume per unit time, Q, of fluid leaking through 
the seals as a function of the applied normal (or squeezing) 
pressure Pq. Results are shown for self afflne fractal surfaces 
with the Hurst exponent H = 0.8 (or fractal dimension Df — 
2.2), and for surfaces with the root- mean-square roughness 
(rms) 1, 2,4 and 6 /im. The fluid pressure difference between 
the two sides is AP = 0.01 MPa and the fluid viscosity /i — 
10"^ Ns/m^ (water). 

of atoms or dynamical variables even for small systems. 
Therefore, we developed a multiscale molecular dynam- 
ics approach to study contact mechanics, to avoid this 
trouble[l2|. The lateral size of the system is L = 1040 A. 
Lx = N^a and Ly = Nya, where a = 2.6 A is the lattice 
space of the block, = Ny ^ 400 for the block. The 
elastic modulus and Poisson ratio are E — 77.2 GPa and 
1/ = 0.42. The lattice space of the substrate is 6 « a/<j>, 
where (/> = [l + a/S) /2 is the golden mean, in order to 
achieve (nearly) incommensurate structures at the inter- 
face. 

For self-affine fractal surfaces, the power spectrum has 
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FIG. 9: The lateral size Ac ~ ^{(c) of the critical constriction 
of the percolation channel, as a function of the applied normal 
(or squeezing) pressure Pq. Results are shown for self afflne 
fractal surfaces with the root-mean-square roughness (rms) 
2 fim and for the Hurst exponent H = 0.9, 0.8, 0.7 and 0.6. 
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FIG. 10: The interfacial separation Uc = 'Ui(Cc) a,t the critical 
constriction of the percolation channel, as a function of the ap- 
plied normal (or squeezing) pressure Po- Results are shown for 
self afflne fractal surfaces with the root-mean-square rough- 
ness (rms) 2 ^m and for the Hurst exponents H = 0.9, 0.8, 
0.7 and 0.6. 

power-law behavior C{q) ~ g-2(if+i)^ where the Hurst 
exponent H is related to fractal dimension Df of the sur- 
face via H — 3 — Df. For real surfaces this relation holds 
only for a finite wave vector region qi < q < qo, where 
qi = 2n/b, is roll-off wave- vector go = Sg^ (see Fig. 
I13p . The randomly rough surfaces have been generated 
as described in Ref. [3, [l2| , which have root- mean-square 
roughness hrms = 10 A and fractal dimension Df — 2.2. 
The roll-off wave- vector go = 3gL, where — 2tt/L and 
L = 1040 A. In this section we define the magnification 

The atoms at the interface between block and substrate 
interact with repulsive potential U{r) = e(ro/r)^'^, where 
r is the distance between a pair of atoms, ro = 3.28 A 
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FIG. 11: The volume per unit time, Q, of fluid leaking 
through the seal as a function of the applied normal (or 
squeezing) pressure Po- Results are shown for self afiine 
fractal surfaces with the root-mean-square roughness (rms) 
2 /im and for the Hurst exponent H = 0.9, 0.8, 0.7 and 
0.6. The fluid pressure difference between the two sides is 
AP = 0.01 MPa and the fluid viscosity /i = 10"^ Ns/m^ 
(water). 



qi 3qL qi 
logq 

FIG. 13; Surface roughness power spectrum of a surface which 
is self-affine fractal for q\ > q > 3qL- The slope logC — 
logq relation for q > 3qL determines the fractal exponent of 
the surface. The lateral size L of the surface determines the 
smallest wave-vector qL ~ 2tt/L. 
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(b) !,=3, A/AO^O.778 
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FIG. 12: The interfacial separation itc — iti(Cc) and the rms- 
roughness /irms(Cc) in the critical constriction of the percola- 
tion channel, as a function of the applied normal (or squeez- 
ing) pressure Pq. Results are shown for self a afiine fractal 
surfaces with the Hurst exponent H = 0.8 (or fractal dimen- 
sion D{ = 2.2), and with the root-mean-square roughness 
(rms) 6 fim. 

and e = 74.4 meV. In molecular dynamics simulations 
there is no unique definition of contact (see Ref. [l^). 
Here we use the critical distance dc to define contact. 
If the separation between two atoms is smaller than dc, 
it has been denoted as contact, otherwise non-contact. 
Here 4 = 4.36 A. 

Fig. [14] shows the block-substrate contact region at 
different magnifications C — 1,3,6,9,12,648. Note that 
when the magnification is increased from 9 to 12, the 
non-contact region percolates. The percolation occurs 
when the normalized projected contact area A/Aq « 0.4, 




FIG. 14: The contact regions at different magnifications 
= 1,3,6,9,12,648, are shown in (a)-(f) respectively. The 
pressure is p « 4.1 GPa. When the magnification is increased 
from 9 to 12, the non-contact region percolate. 



in good agreement with percolation theory [11|. 

5. Improved analytical description 

In Sec. 2 we assumed that all the fluid flow occurs 
through a single constriction which we refer to as the 
critical constriction. In reality, fluid flow will also occur 
in other flow channels even if they have more narrow con- 
strictions. In this section we will assume that there are 
a finite concentration of critical or near critical constric- 
tions, which correspond to all constrictions which appear 
when the magnification changes in some narrow inter- 
val around the critical value Cc, e.g., in such a way that 
A{Q/Aq changes by, say, ±0.03. Since we are very close 
to the percolation threshold, we will assume that the size 
of all the (nearly critical) constrictions remains the same. 
In a more accurate treatment one would instead intro- 
duce a distribution of sizes of constrictions. In Fig. [TSTal 
the dots correspond to the critical or near critical con- 
strictions along percolation channels (solid lines). One 
expects the (nearly critical) constrictions to be nearly 
randomly distributed in the apparent contact area, and 
that the channels, of which they are part, to have all pos- 
sible directions as indicated by the lines in Fig. fTsT a). 
Here we will consider a simplified version of (a) where 
the (nearly critical) constrictions form a more ordered 
arrangement as in Fig. llSf b). In reality, the dots and 
the lines should be (nearly) randomly distributed as in 
llSf a) , but this is likely to have only minor effects on what 
follows. 

On the average fluid will only flow in the x-direction. 
Thus, in a first approximation one may assume that 
no fluid flows along the (transverse) channels pointing 
(mainly) in the j/-direction in Fig. [TSFb). Let a be the 
(average) distance between two nearby critical constric- 
tions. Thus we expect n = L^/a constrictions along a 
percolation channel (in the figure we have n = 3). If Qi 
denotes the fluid volume per unit time flowing along one 
percolation channel, then we must have 

Qi = M{P,-Pi) = M{Pi-P2) = ... = MiPr,-P\^. (5) 

From (5) we get 

Ql = — (^a-Pb) 

n 

As expected, the amount of fluid flowing in the channel is 
reduced when the number of constriction increases. How- 
ever, there will be roughly Ly/a percolation channels so 
the total fluid flow will be 

Q ^ bLM{P, - Pb) (6) 

which is identical to the result obtained in Sec. 2. This 
analysis is very rough, and a more detailed analysis will 
result in some modifications of the leak-rate, but (6) 
should be very useful as a first rough estimate of the 



(a) 




constriction 

FIG. 15: The solid lines denote non-contact channels and 
the dots critical or near critical constrictions. In reality the 
constrictions and channels are nearly randomly distributed as 
in (a) [see also Fig. I14f e)] but in the model calculation we 
use the more ordered structure shown in (b). 
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FIG. 16: Nominal contact pressure distribution Po{x) (curve 
bounding the dotted area) for (a) O-ring seal and (b) lip seal. 
In (a) the curve denoted by "approx" is an approximation to 
the continues Hertz-like curve. 

leak-rate. Note that the present treatment will result in 
a more gradual changes in the liquid pressure in the ap- 
parent contact region, from the initial high pressure value 
Pa (entrance side) to the low pressure value Pb (exit side). 

6. Comparison with experiment 

We have not found any results in the literature about 
leak-rates of seals for well characterized systems. How- 
ever, we have found some results which are in qualitative 
agreement with our theory. For example, leak-rates ob- 
served for both rubber and steel seals tend to decrease 
very fast (roughly exponentially) with increasing contact 
force. Thus, in Ref. the leak-rate for a rubber seal 
decreased by 6 order of magnitude as the load increased 
by a factor of 10. A similar sharp drop in the leak-rate 



with increasing contact force has been observed for seals 
made from steel [31 ■ However, in the latter case some 
plastic deformation is likely to occur in the contact re- 
gion. In both cases the nominal pressure may change 
less than the change in the load, due to an increase in 
the nominal contact area with increasing load. A de- 
tailed analysis of the experimental data is not possible as 
the surface topography was not studied in detail. In Sec. 
9 we suggest a simple experiment which can be used to 
test the theory. 

The present theory implys that most of the fluid leak- 
age occurs through the critical or nearly critical constric- 
tions in the percolating channels at the interface between 
the two solids. Since the constrictions are very small they 
can easily be clogged up by dirt particles in the fluid. 
This results in leak-rates which decreases with increas- 
ing time as the microscopic gaps get clog ged up. This 
has recently been observed for metal seals IJ]. In fact. 



by using specially prepared fluids with immersed parti- 
cles having a narrow distribution of particle diameters, it 
should be possible to determine (approximately) the size 
(or rather the height) of the critical constriction. 

7. Comment on the role of non-uniform pres- 
sure, rubber viscoelasticity and adhesion 

In the study above we have assumed that the normal 
(squeezing) pressure is constant in the nominal rubber- 
countersurface contact region. In reality, this is (almost) 
never the case. Thus, in rubber 0-ring applications one 
expect a pressure distribution which is Hcrzian-like, as 
indicated in Fig. [TBTa). In (dynamical) rubber seals 
for linear reciprocal motion, the pressure distribution is 
asymmetric, with a much steeper increase in the pres- 
sure when going from the high-pressure (Pa) fluid side 
towards the center of the seal, as compared going from 
the low-pressure (Pb) fluid side toward the center of the 
seal, see Fig. [T6lb). (The reason for this asymmetry 
does not interest us here.) The theory developed above 
can be applied approximately to these cases too. Thus in 
case (a) (e.g., rubber 0-ring seals) one may approximate 
the actual Herzian-like pressure profile with a sum of step 
functions as indicated in Fig. [TBJa). Since the seal-action 
is so strongly dependent on the squeezing pressure (see 
Figs. [8l and [TT|) . it is enough to include the central step 
region (with width L) in the analysis. Since there is no 
unique way to determine the width L there will be some 
(small) uncertainty in the analysis, but this is not im- 
portant in most practical cases. Similarly, for the lip-seal 
[Fig. I16r b)] during stationary condition it is enough to 
include the region (width L) where the normal pressure 
is maximal. 

In the study above we have assumed that the rub- 
ber behaves as a purely elastic solid. In reality, rubber- 
materials are viscoelastic. One consequence of this is 
stress relaxation. For example, after a rubber 0-ring has 
been deformed to fit into the "cavity" where it is placed, 



the stress exerted on the solid walls will decrease with 
increasing time. Since rubber materials have very wide 
distribution of relaxation times, the stress can continue 
to decrease even one year after installation. Thus, after 
very long time the pressure in the rubber-countersurface 
contact region may be so low that the seal fails (note: we 
found earlier that the leak-rate depends extremely sensi- 
tively on the normal pressure). Stress relaxation can be 
easily taken into account in the analysis above by using 
the relaxation modulus E{t) (where t is time) measured 
in the laboratory using standard methods. 

Finally, let us comment on the role of adhesion in rub- 
ber seals. We first note that if the fluid is an oil, the ef- 
fective adhesion between the rubber and the hard coun- 
tersurface may vanish or nearly vanish, as observed in 
some experiments [isjl . If the fluid is not an oil (e.g., wa- 
ter) some effective adhesive interaction may remain. In 
particular, if the fluid is a gas then the the effective ad- 
hesion may be similar to that in the normal atmosphere. 
However, even in this case the adhesive interaction be- 
tween the solids may have a negligible influence on the 
leak-rate. The reason for this is that adhesion operates 
mainly at very short length scales, corresponding to high 
magnification > (^ad, while the leak-rate is determined 
mainly by the contact mechanics at the point where the 
first percolation channel appears, corresponding to the 
magnification (^c- If Cc << Cad the adhesive interaction 
will have a negligible influence on the leak-rate. We now 
illustrate this with a numerical calculation using the the- 
ory of Ref. 01 ■ 

Fig. [TTI shows the relative contact area A{()/Ao as a 
function of the logarithm of the magnification (. Note 
that at the magnification Cc, where the non-contact area 
first percolates, the adhesional interaction has no influ- 
ence on the contact area. The adhesional interaction will 
manifest itself only for C > Cad, where the adhesional in- 
teraction increases the contact area as compared to the 
case without the adhesional interaction included. The 
results in Fig. [T7| is for rubber block in contact with a 
hard solid with a self-affine fractal surface with the root- 
mean-square roughness h^ms = 6 /xm, the Hurst exponent 
H — 0.8, and for the squeezing pressure Pq = 0.2 MPa. 

Finally we note that if there is very little fluid at the 
interface strong capillary adhesion may occur between 
the surfaces. This is known to be of great importance in, 
e.g., the context of rubber wiper blades. This topic has 

flEl 



been discussed in detail in Ref. 



8. Dynamical seals 

The theory presented above is for static seals. Here 
we give some comments related to dynamical seals. We 
will estimate the leak-rate for linear reciprocal seals at 
very low sliding velocity. We assume that the roughness 
occurs mainly on the rubber surface and we treat the 
hard countersurface as perfectly flat. Thus, as the rub- 
ber slides along the countersurface the contact mechan- 
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FIG. 17: The relative contact area A{Q/Ao as a function 
of the logarithm of the magnification (. Note that at the 
magnification ^c, where the non-contact area first percolate, 
the adhesional interaction has no infiuence on the contact 
area. The adhesional interaction will manifest itself only for 
C > Cad, where the adhesional interaction increases the con- 
tact area as compared to the case without the adhesional in- 
teraction included. For rubber block in contact with a hard 
solid with a self-affine fractal surface with the root-mean- 
square roughness /irms = 6 fj,m, the Hurst exponent H = 0.8. 
The squeezing pressure Pq = 0.2 MPa and, for the curve "with 
adhesion", with the interfacial binding energy per unit area 
A7 = 0.05 J/m^ 



ics does not change, e.g., the percolation channel will be 
time-independent in the reference frame moving with the 
rubber. We consider the system in the reference frame 
where the rubber is stationary while the hard countersur- 
face moves from left to right with the velocity vq. The 
rubber is assumed to be below the countersurface, see 
Fig. 1181 The high pressure fluid region (pressure Pa) oc- 
cupys a; < while the low pressure region (pressure Pb) 
occupys X > L. 

We assume a Newtonian fluid and stationary and lam- 
inar flow. The basic equations for the fluid flow are 



Vp — 77V V, 



V-v = 0, 



where p(x) and v(x) are the fluid pressure and the fluid 
flow velocity, respectively. We now consider the fluid 
flows in the percolation channel. Let s be the length- 
coordinate along the percolation channel. Since in gen- 
eral the A(s) >> 'w(s), where A(s) is the width and u{s) 
the height of the channel at the point s along the channel, 
we can write the velocity as v(x) = sv{s, z) where 

/ \ ^ dp ^ ^ z 

v{s, z) w — — z{z - u{s)) + VqX ■ s—— 
Irj as u[s) 

The volume flow per unit time through any cross section 
of the channel is assumed to be the same, and equal to 
Q which gives 

nu(s) 

Q = A(s) / dz v{s) = 
Jo 



or 



w ^ , u^js) dp . Ms) 



1277 ds 
dp 677 

-VqX ■ S ■ 



ds u'^{s) X{s)u^{s) 
Integrating this equation gives 



P{Q = Pa + / ds 
JQ 



6?7 

u^(s) 



VqX ■ S - 



12?7Q 

X{s)u^s) 



(7) 



where l^, is the length of the percolation channel. Let 
P(u) be the probability to flnd the surfaces separated by 
the hight u along the percolation channel. Note that A(s) 
also can (at least locally) be considered as a function of 
u, which we denote by A(u) for simplicity. Thus we can 
write (7) as 



p{l 



Pa+QL^rivo / du 



Piu) 



du 



Piu) 



(8) 

where La is the length of the percolation path projected 
on the X-axis. In Ref. [§] we have shown how it is possi- 
ble to calculate the distribution P„ of heights u between 
two surfaces in elastic contact. We will now assume that 
(note: u > Uc) 

P{u) 



I^du'P^, 
Let us write (8) as 

p{Q =P'^^P^ + B,vo - Ca 



where 



and 



Pa = 6Pa'7 



du 



P{u) 



^du 



(9) 



(10) 



(11) 



Similarly, one getsfisj 

p{lb) = P^ = Pb - Pb«o - CbQ (12) 
Thus in this case (1) takes the form 
Q = M{Pi - P^) = 

M(Pa - Pb) + A/(Pa + Bb)wo - M(a - Cb)g 



or 



AP + (Pa + B^)vo 
1 - M(a - Cb) 



(13) 



z 



The factor M{Ca — Cb) in the denominator in this ex- 
pression is independent of vq and we will assume that 
it is negligible compared to unity, and neglect it. One 
interesting application of (13) is to wiper blade. Here 
AP = so that (13) takes the form 



Q = MiB, + B,,)vo 



(14) 



Substituting (2) and (10) (and a similar expression for 
Bb) in (14) gives 



du 



P{u) 



(15) 



where we have included the extra factor Ly/L^ to take 
into account the number of square seal units. During 
the time t the leak-volume is Qt. We define the average 
thickness d of the leak-film as d = Qt / (LyVot) . From (15) 
we get 



d — f3uc 



^-1 



du 



,P{u) 



(16) 



(17) 



We have calculated the integral / in /3 for some typical 
casesfl^, and found that / w 0.1 — 0.2 so we expect 
0.1. 

In the treatment above we have assumed that the con- 
tact between the rubber and the hard countersurface does 
not depend on the pressure in the fluid, which is a good 
approximation as long as the fluid pressure p(x) << Pq- 
However, as the sliding velocity increases, the fluid pres- 
sure in some regions at the interface will increase which 
will tend to increase the separation between the two sur- 
faces. At very high sliding velocity, hydrodynamic lubri- 
cation will prevail and the surfaces are completely sepa- 
rated by a thin fluid fllm. However, even at much lower 
sliding velocity the hydrodynamic pressure buildup may 
strongly increase the leak rate. In particular, the pres- 
sure at the critical constriction will tend to increase the 
separation between the surfaces and hence increase the 
leak-rate. We will not study this effect here but just es- 
timate when this effect becomes important. Let Pc be 
the pressure at the critical constriction. If pc acts over 
the area it will locally increase the separation between 
the surfaces by the amount 20] Au XcPc/E. Thus, the 
pressure at the critical constriction must be much smaller 
than Euc/Xc in order for the pressure-induced effect to 
be negligible. Note that 



Pc 



GLrivo 



du 



Piu) 



u^ 



so the present study is limited to sliding velocities 

-1 



vo « 



GLrjXc 



du u, 



Eul 
Lr]Xc 



(17) 




FIG. 18: A rubber block with a rough surface in contact 
with a hard smooth countersurface (upper block) which moves 
relative to the rubber block with the velocity vq. 



where we have used that the integral typically is of or- 
der ~ 0.15. Thus, for example, if in a wiper blade 
application [2l|. after some use the rubber blades typi- 
cally develop (because of wear) a surface roughness with 
a rms amplitude of several micrometer. If the rms rough- 
ness is 2 fim (and the Hurst exponent H = 0.8), the 
nominal pressure ~ 0.2 MPa, and if we assume that 
E ^ 10 MPa we get from Fig. [7]and[5]uc ~ 0.1 /im and 
Ac ~ 10 /im. If L w 0.1 mm and (for water) w lO^'^ Pas 
we get that the slip velocity must be at most ~ 1 cm/s 
in order for (17) to be valid. According to (16) the (av- 
erage) fllm thickness of the water layer would be of order 
0.01 nm. 

9. A new experiment 

Very few studies of leak rates of seals with well char- 
acterized surfaces have been published. Here we would 
like to suggest a very simple experiment which could be 
used to test the theory presented in Sec. 2. In Fig. [TOlwe 
show a set-up for measuring the leak-rate of seals. A glass 
(or PMMA) cylinder with a rubber ring (with rectangu- 
lar cross-section) glued to one end is squeezed against a 
hard substrate with well-deflned surface roughness. The 
cylinder is flUed with a fluid, e.g., water, and the leak- 
rate of the fluid at the rubber-countersurface is detected 
by the change in the height of the fluid in the cylinder. In 
this case the pressure difference AP = P^ — P^ = pgH, 
where g is the gravitation constant, p the fluid density 
and H the height of the fluid column. With « 1 m 
we get typically AP fa 0.01 MPa. With the diameter 
of the glass cylinder of order a few cm, the condition 
Pq >> AP (which is necessary in order to be able to 
neglect the influence on the contact mechanics from the 
fluid pressure at the rubber-countersurface) is satisfled 
already for loads (at the upper surface of the cylinder) of 
order kg. 

10. Summary and conclusion 

Seals are extremely useful devices to prevent fluid leak- 
age. However, the exact mechanism of roughness induced 



improved and extended in various ways. 
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FIG. 19: A simple experimental set-up for measuring the leak- 
rate of seals. A glass (or PMMA) cylinder with a rubber 
ring glued to one end is squeezed against a hard substrate 
with well-defined surface roughness. The cylinder is filled 
with a fiuid, e.g., water, and the leak-rate of the fluid at the 
rubber-countersurface is detected by the change in the height 
of the fluid in the cylinder. In this case the pressure difference 
AP = — Ph = pgH, where g is the gravitation constant, p 
the fluid density and H the height of the fluid column. With 
?a 1 m we get typically AP ^ 0.01 MPa. 



leakage is not well understood. We have presented a the- 
ory of the leak-rate of seals, which is based on percolation 
theory and a recently developed contact mechanics the- 
ory. We have studied both static and dynamics seals. We 
have presented numerical results for the leak-rate Q, and 
for the lateral size Ac and the height Uc of the critical 
constriction. We assumed self affine fractal surfaces and 
presented results for how Q, Ac and Uc depend on the 
root-mean-square roughness amplitude and the fractal 
dimension D f = 3 — H (where H is the Hurst exponent), 
and on the pressure Pq with which the rubber is squeezed 
against the rough countersurface. 

We have also presented molecular dynamics results 
which show that when two elastic solids with randomly 
rough surfaces are squeezed together, as a function of in- 
creasing magnification or decreasing squeezing pressure, 
a non-contact channel will percolate when relative pro- 
jected contact area, A/Aq, is of order 0.4, in accordance 
with percolation theory. Finally, we have suggested a 
simple experiment which can be used to test the theory. 

The theory we have presented in this paper is very 
rough, but we believe that it captures the most impor- 
tant physics, and that the presented approach can be 



Acknowledgments: We thank Ed Widder (Federal 
Mogul Scaling Systems) and Matthias Schmidt (IFAS, 
RWTH Aachen University) for useful correspondence re- 
lated to seals. We also thank A. Koenen (Valeo Systeme 
d'Essuyage) for very interesting information related to 
the tribological properties of the contact between rubber 
and glass in the context of rubber wiper blades. 



[4; 

[5 
[G 
[7 

[s: 
[9: 

[10 

[11 

[12 

[13 
[14 

[15 

[16 

[17 

[18 



[19] 



R. Flitney, Seals and sealing handbook (Elsevier, 2007). 
M. Mofldi, B. Prakash, B.N.J. Persson and O. Albohr, J. 
Phys.: Condens. Matter 20, 085223 (2008). 
See, e.g., B.N.J. Persson, O. Albohr, U. Tartaglino, A.I. 
Volokitin and E. Tosatti, J. Phys. Condens. Matter 17, 
Rl (2005). 

J. Chem. Phys. 115, 3840 (2001). 
Phys. Rev. Lett. 99, 125502 (2007). 
Surf. Science Reports 61, 201 (2006). 
Eur. Phys. J. E8, 385 (2002). 
F. Bucher and B. Chiaia, Phys. Rev. 



Persson, 
Persson, 
Persson, 
Persson, 
Persson, 



B.N.J. 
B.N.J. 
B.N.J. 
B.N.J. 

B. N.J. 
B65, 184106 (2002). 

C. Yang and B.N.J. Persson, J. Phys.: Condens. Matter, 
in press 

B. N.J. Persson, O. Albohr, C. Creton and V. Peveri, J. 
Chem. Phys. 120, 8779 (2004) 

D. Stauffer and A. Aharony, An Introduction to Percola- 
tion Theory, CRC Press (1991). 

C. Yang, U. Tartaglino and B.N.J. Persson, Eur. Phys. J 
E19, 47 (2006). 

E. Widder, ASME B46 Seminar, AprU 15, 2004. 

M. Schmidt, H. Murrenhoff, H. Lohrberg and F-J. 
Korber, to be published. 

B. Zappone, K.J. Rosenberg and J. Israelachvili, Trib. 
Lett. 26, 191 (2007). 

B.N.J. Persson, submitted to Journal of Physics: Con- 
dens. Matter. 

F. Deleau, D. Mazuyer and A. Koenen, Tribology Inter- 
national (in press). 

We assume that no cavitation occur at the exit of the 
critical constriction. Cavitation may occur if the local 
pressure at the exit of the critical constriction is nega- 
tive, but in this paper we assume so low sliding velocity 
vo that this is not the case. If cavitation occur p{lh) ~ 0, 
or, more accurately, p{lh) is equal or close to the vapor 
pressure of the fluid or of gases dissolved in the fluid. 
We note that the transition from boundary or mixed lu- 
brication to hydrodynamic lubrication probably involves 
cavitation at the exit of many narrow constrictions, since 
otherwise the total load supported by the fluid film will 
be very small. It would be very interesting to study this 
problem theoretically, since the transition from bound- 
ary lubrication to hydrodynamic lubrication is not well 
understood. 

We did two calculations: For a self afline fractal surface 
with the rms roughness of ftrms = 2 fim and the squeezing 
pressure Po = 0.18 MPa we got B « 0.11 (and Uc ~ 
0.1 nm and Ac ~ 16 /xm) and for /irms ~ 4 /im and 



P() = 0.2 MPa we got B « 0.21 (and Wc « 1 Mm and 
Ac ~ 109 Mm). 

] The fluid pressure along the pereohition ehannel in- 
creases towards the critical constriction. Thus, the elas- 
tic deformation of the rubber at the critical constric- 
tion is determined not just by the pressure at the con- 
striction but also by the pressure acting on the rubber 
along the percolation path. If one assume a strait per- 



colation path and that the pressure along the percola- 
tion path p{x) = pcx/L then one can easily show that 
Au ~ ln(L/Ac)AcPc/-B but the additional logarithmic 
factor is never very large and does not change our quali- 
tative conclusion. 
[21] A. Koenen and A. Sanon, Tribology International 40, 
1484 (2007). 



